5d-5f Electric- multipole Transitions in Uranium Dioxide Probed by Non-resonant 

Inelastic X-ray Scattering 
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Non-resonant inelastic x ray scattering (NIXS) experiments have been performed to probe the 5d- 
5/ electronic transitions at the uranium 04,5 absorption edges in uranium dioxide. For small values 
of the scattering vector q, the spectra are dominated by dipole-allowed transitions encapsulated 
within the giant resonance, whereas for higher values of q the multipolar transitions of rank 3 and 
5 give rise to strong and well-defined multiplet structure in the pre-edge region. The origin of the 
observed non-dipole multiplet structures is explained on the basis of many-electron atomic spectral 
calculations. The results obtained demonstrate the high potential of NIXS as a bulk-sensitive 
technique for the characterization of the electronic properties of actinide materials. 

PACS numbers: 78.70. Ck, 71.70-d 



Non-resonant inelastic x ray scattering (NIXS) is 
emerging as a powerful bulk-sensitive probe of the dy- 
namical electron-density response of solids and as a new 
tool for the investigation of the electronic structure in 
strongly-correlated systems . The sharp non-dipolar 
pre-thrcshold features that dominate the NIXS energy- 
loss spectrum at large momentum transfer Q can be used 
to study atomic-environment, valence, and hybridization 
effects [6| and give information on states with different 
symmetries to those of electric-dipole spectra obtained 
by x ray absorption spectroscopy (XAS). No NIXS ex- 
periments have been performed so far on actinide solids, 
and the main aim of this work has been to explore the 
potential of NIXS for studying the nature of the 5/ elec- 
tron shell through the observation of multipole transi- 
tions from core 5d to valence 5/ states. The actinide 
series is still only modestly understood. Several funda- 
mental questions remain unanswered, including the num- 
ber of electrons in valence states, the angular momentum 
coupling, and the character of the bonding. These factors 
have a huge influence in driving the peculiar magnetic 
and electronic properties of actinide metals Q, alloys, 
and compounds [H, and the development of novel spec- 
troscopic techniques sensitive to the ground-state prop- 
erties of these systems is of high importance. 

Compared to resonant inelastic x ray scattering (RLXS) 
and electron energy loss spectroscopy (EELS), the NIXS 
technique is penalized by a lower intensity, but it has the 
advantage that it permits a straightforward and quanti- 
tative interpretation of the spectra. Element specificity 
is provided by the energy positions of the core to va- 
lence transitions. Moreover, even when interrogating 
low-energy edges, NIXS remains a bulk-sensitive probe 
since the sampling depth is largely determined by the 
incident photon energy. The high penetration depth of 



hard X-rays is also advantageous if the sample must be 
inside an outer container, e.g., if thermodynamical pa- 
rameters such as temperature and pressure must be var- 
ied, or when the sample must be encapsulated for safety 
reasons. 

Here we report the results of NIXS studies on UO2, a 
widely investigated semiconducting compound exhibit- 
ing long-range order of magnctic-dipole and electric- 
quadrupole moments at low temperature @-[l0j]. The 
U 4+ ions have a 5/ 2 configuration with a 3 #4 ground 
state, with localized 5/ states in the 5 eV gap between the 
O 2p band and empty U 6d band, with a 5/ 2 — > 5/ 1 6sd 
excitation energy of 2-3 eV [IlJ . 

The experiment has been performed using the ID16 
inverse-geometry, multiple-analyzer-crystal spectrome- 
ter at the European Synchrotron Radiation Facility in 
Grenoble, France 12J. The beam generated by three con- 



secutive undulators was monochromatized by a Si(lll) 
double-crystal, and horizontally focussed by a Rh-coated 
mirror. The beam size at the sample position was 0.3x 1.3 
mm 2 (horizontal x vertical). A set of nine spherically- 
bent Si(660) analyzer crystals with 1 m bending radius, 
horizontal scattering geometry and vertical Rowland cir- 
cles, provided a bandwidth of 1.3 eV at a final photon 
energy Ef of 9.689 kcV, and an intensity of 7xl0 13 pho- 
tons/s for a 25 fj, rad vertical divergence of the undulator 
radiation. The Bragg angle of the analyzers was fixed 
at 88.5°, and the geometry chosen in order to measure 
five different momentum-transfer values simultaneously, 
using a position sensitive detector based on a 256x256 
photon counting pixel array. For the measurements we 
used a UO2 single crystal with an external surface per- 
pendicular to a [111] direction. A He-flow cryo-cooler 
gave us the possibility to cool the sample to a base tem- 
perature of 10 K. Data have been collected in symmet- 
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ric reflection geometry by scanning the incident energy 
Ei = Ef + hiu at fixed final energy Ef, covering the ftw 
energy-loss interval of the uranium 5 (5d 5 / 2 — > 5/) and 
O4 (5d 3 /2 5/) absorption edges. Typical counting time 
was six hours per spectrum. Cross sections are particu- 
larly large between shells of the same principal quantum 
number, such as for the 5d — > 5/ transition. We also per- 
formed measurements across the uranium (4c? — > 5/) 
absorption edges, but we failed to observe NIXS signals 
at these edges due to the much lower cross section. 

The double differential cross section for NIXS events 
can be expressed within the first Born approximation and 
is proportional to a dynamic structure factor 5(q, uj) re- 
lated to the complex dielectric-response function e(q, oS) 
by the fluctuation dissipation theorem [l3, 14|, S(q, lo) oc 
q 2 Im[e(q, w) -1 ]. The non-resonance radiation-matter in- 
teraction is dominated by a term proportional to the 
square of the vector potential A and is given by 
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kf is the scattering vector, eij and 
k.j j arc respectively the polarization and wave vectors of 
the incident and scattered radiation. By expanding the 
transition operator e lq r in terms of spherical harmonics 
Eq. (TTJ) can be written as @, @] 
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where the factors Dk, m describe the angular dependence 
of the signal and jk{qr) are spherical Bessel functions 
of rank k. The multipole moments k for the I — > (' 
shell transition are limited by the triangular condition, 
<k< l + l', and the parity rule, l + l' + k = even. 
Thus for d / transitions, only k = 1 (dipolc), k = 3 
(octupolc), and k = 5 (triakontadipolc) transitions are 
allowed. The dipole transition dominates the response of 
the system for low values of q, whereas at higher q values 
the higher multipole transitions become important. The 
different 2 fc -pole transitions excite the system to different 
subsets of different symmetries and energies within the 
final state configuration. 

The NIXS spectra obtained for UO2 at room temper- 
ature are shown in Fig. [1] The spectrum measured for 
q = 9.6 A -1 over the whole energy- loss range investi- 
gated is shown in the inset. The NIXS signal has been 
obtained by fitting the broad Compton peak to a polyno- 
mial function and by subtracting it from the total mea- 
sured intensity. The data at different q values have been 
normalized to the peak intensity of the feature at about 
104 eV. The dipole spectrum observed at q = 2.81 A -1 
shows the typical, ill-defined giant-resonance shape to- 
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FIG. 1: (Color online). NIXS spectra measured for UO2 at 
the uranium 04,5 edges with a fixed final energy of 9.689 
keV, at different values of the scattering vector q. The inset 
shows the spectrum measured in an extended energy range, 
including the elastic and Compton peaks. The NIXS signal 
has been normalized to the peak at 104 eV in the high-q 
spectra. 



gether with a prepcak around 99 eV, in good agree- 
ment with earlier XAS [l6[ and EELS 17, 1^ measure- 
ments. Increasing q above 9 A -1 , dipole transitions be- 
come negligible while higher multipole transitions appear 
at lower energies. An asymmetric peak with a maximum 
at 103.6 eV and a clearly-resolved double-peak feature 
with maxima at 94.9 and 97.3 eV characterize the spec- 
trum at q = 9.88 A~F NIXS spectra for UO2 have also 
been measured at 10 K, in the ordered antifcrromagnetic- 
antifcrroelectric quadrupolc phase. A comparison with 
the room temperature spectra does not reveal differences 
beyond the statistical error. 

Many-electron atomic spectral calculations enable us 
to identify the origin of the observed multipole transi- 
tions. The initial- and final-state wave functions have 
been calculated in intermediate coupling using Cowan's 
atomic Hartree-Fock (HF) code with relativistic correc- 
tions [ljj], which gives the Slater parameters, Fjj = 
9.711, Ff f = 6.364, F 6 ff = 4.677, Fj f = 10.652, Ff 



ff 

6.850, G\ f = 12.555, G 3 df = 7.768, G b d} = 5.544, and spin- 
orbit parameters, Cs/ — 0.274, ^5rf = 3.199 (all values 
in eV). Angle-integrated electric 2 fc -pole transitions have 
been calculated from initial state U 5/ 2 to the final-state 
levels of the U 5d 9 5/ 3 configuration in spherical sym- 
metry, implicitly using the selection rules and symmetry 
constraints. The calculated line spectra have been broad- 
ened by a Lorentzian of T = 0.2 eV and a Gaussian of 
a = 0.4 eV. The results are shown in Fig. [21 where the 
Slater parameters were scaled to 70% of the HF values to 
account for screening effects. The decay processes leading 
to the giant resonance were not included in the calcula- 
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FIG. 2: (Color online). Top panel: 

tra for UO2 measured at 1.68 A _i (red dashed line) and 
9.74 A -1 (blue solid line) compared to results of atomic mul- 
tiplet calculations for the U 5/ 2 configuration (lower panel). 
Calculated contributions to the cross section from dipole (k = 
1), octupole (k = 3), and triakontadipole (k = 5) transitions 
are shown separately. 



tion. This affects the broadening of dipole spectrum, as 
will be discussed below. 

The cross sections of the three different channels 
change gradually as a function of q. At q = 2.81 A -1 
only dipole transitions play a role, but with increasing q 
the higher multipole transitions become prominent. The 
ratio of the octupole to triakontadipole cross sections is 
~0.6 and ~0.85 at q = 9.15 and 9.88 A -1 , respectively, 
with negligible dipole cross section. 

Fig. Q] shows that for q = 2.81 A -1 the dipole-allowed 
peaks lie within the energy region of the giant resonance. 
The prepeak at ~99 eV has previously been attributed 
to the finite spin-orbit interaction 17J, LL8[ • Going from 
q = 9.15 to 9.88 A -1 , the experimental data in Fig. [1] 



show changes that are in agreement with the calculated 
spectra, noting that with increasing q the k =5 contribu- 
tion increases with respect to the k = 3. The intensity of 
the low-energy peak in the doublet structure (~95 and 
97 eV) increases relative to that of the high-energy peak 
with increasing q. Concurrently, the asymmetric peak 
at ~104 eV increases in intensity and also shifts slightly 
to lower energy, confirming that it contains two different 
k contributions. Thus we can conclude that there is a 
good overall agreement between experimental and calcu- 
lated spectra. In order to explain the differences between 
the multipole spectra and to understand the pronounced 
multiplet structure we need to go into the details of the 
selection rules and the interplay between electrostatic 
and spin-orbit interactions. In NIXS, the final-state J' 



and L' values are restricted by|J — k\ < J' < J + k and 
\L — k\ < V + AS < L + k, respectively, where AS = 
in LS coupling and AS = ±1 in intermediate coupling. 
Thus for initial state 5/ 2 3 i?4, dipole transitions are al- 
lowed to final states 3 G, 3 H, 3 I with J' = 3, 4, 5, where 
in intermediate coupling quintet-spin states can mix in 
by 5d core spin-orbit interaction. The selection rules have 
implications for the final states that can be reached. For 
the transitions /" — > d 9 f n+1 with less than half-filled 
shell, the maximum spin in the initial state is one less 
than in the final state. Since electric-multipole transi- 
tions obey the selection rule AS = 0, the high-spin final 
state is a ghost state, which cannot be reached from the 
initial state without spin-orbit interaction (20j . A simi- 
lar situation arises for L' and J' . While the / 2 ground 
state has L = 5 and J = 4, the final state d 9 f 3 contains 
values up to L' = 8 and J' = 7. High L' and J' values 
cannot be reached by dipole transitions, only by higher 
multipole transitions. 

Next, we must consider the energy distribution of the 
allowed final states. The electrostatic interactions give 
an energy spread over the different L'S' states of about 
30 eV (see Fig. [3]). In general, states with higher 5" 
and L' values have lower energies; however, it must be 
stressed that the second Hund's rule is not strictly valid 
in the final state, so that states with intermediate L' 
value can have a lower energy [2(| . Fig. [3][d) shows that 
in LS coupling, where we only take the electrostatic in- 
teractions into account, the dipole-allowed final states 
have much higher energies than the other states. The 
reason is that if the momentum transfer in the excitation 
process is small, the charge cloud in the final state has to 
be oriented in approximately the same way as in the ini- 
tial state, which results in a high electrostatic energy. In 
contrast, higher multipole transitions allow final states 
where the angular distribution of the charge density is 
more evenly spread, thereby reducing the electrostatic 
energy. 

Generally, L'S' states mixed by spin-orbit inter- 
action will have an energy separation of AE = 
^(Adcctrostatic) 2 + ( A spin _ or bit) 2 ■ In the case of the 
dipole-allowed L'S' states with triplet spin, located at 
higher energies, the spin-orbit coupling switches on the 
prepeak at ~99 cV with quintet-spin character 
Fig. (3Ja-b) shows that in the dipole spectrum both 
peak position and prepeak intensity are most sensitive 
to scaling of the electrostatic interactions. On the other 
hand, the higher multipole allowed L'S' states are closely 
grouped together so that they experience a strong spin- 
orbit interaction. Consequently, these states will be 
split into two groups with an energy separation of about 
|C5d ~ 8 eV; one group with core hole 5d 5 / 2 and the other 
with 5d 3 / 2 character (see Fig. [3]). Thus the core hole spin- 
orbit interaction manifests itself prominently by the main 
energy splitting in the higher multipole spectra. As a fi- 
nal point, we discuss the giant resonance appearing as a 
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FIG. 3: (Color online). Calculated multipole spectra with 
k = 1 (black dotted line), k = 3 (blue solid line), and k — 5 
(red dashed line) for the transition U 5/ 2 — > 5d 9 5/ 3 with dif- 
ferent scaling factors g for the Slater integrals and z for the 
spin-orbit interaction, (a) Realistic values: g = 0.7, z — 1 
(intermediate coupling scheme applies), (b) strongly reduced 
electrostatic interactions: g — 0.3, z = 1, (c) only spin-orbit 
interaction (jij-coupling scheme applies), and (d) only elec- 
trostatic interactions (LiS-coupling scheme applies). The ap- 
proximate peak assignment is indicated when meaningful. As 
states are strongly mixed in intermediate coupling peak as- 
signment is less meaningful. 



broad band around 112 eV. This feature appears clearly 
in the measured dipole spectrum, (Fig. [TJ), and is due 
to resonant decay into continuum states. Remarkably, 
no broad structures are observed in the higher multi- 
pole spectra, even though some excited states are close 
to the onset of the giant resonance. A simple explana- 
tion would be that the continuum states lie energetically 
near the dipole-allowed transitions. An alternative ex- 
planation might be that the continuum states have the 
same symmetry as the dipole-allowed states, and hence, 
interfere strongly, but have a different symmetry than 
the higher multipole-allowed states. Furthermore, low 
ener gy s tates with high L' values are often double forbid- 
den 18[, meaning they are not only forbidden in dipole 
excitation, but also do not have decay channels. This 
leads to long lifetimes and, hence, narrow peak widths. 

In conclusion, we have demonstrated the potential 
of the NIXS technique for studying bulk actinide com- 
pounds. The intrinsic intensity limitation of the tech- 
nique can be overcome by using an optimized spectrom- 
eter at a third-generation synchrotron-radiation source. 
The instrumental energy resolution allows us to resolve 
the multiplet structure. For the localized 5/ electrons in 



UO2 a good agreement is obtained with many-electron 
calculations in intermediate coupling. Electric-multipolc 
transitions from the ground state can reach only a lim- 
ited subset of final states, thereby providing a fingerprint 
for the specific ground state. These results are encour- 
aging for perusing further NIXS measurements in a wide 
range of investigations on transuranium materials, with 
a focus on the dual localized-delocalized nature of the 5/ 
electrons and their instability induced by external pres- 
sure. 



[1 

[2 

[3 

K 
[5 
[6 

[7] 
[8] 

[9 
[10 

[11 
[12 

[13 
[14 

[15 

[16 

[17 

[18 

[19 

[20 



A. T. Macrander, P. A. Montano, D. L. Price, V. I. Kush- 
nir, R. C. Blasdell, C. C. Kao, and B. R. Cooper, Phys. 
Rev. B 54, 305 (1996). 

I. G. Gurtubay, J. M. Pitarke, W. Ku, A. G. Eguiluz, 

B. C. Larson, J. Tischler, P. Zschack, and K. D. Finkcl- 
stein, Phys. Rev. B 72, 125117 (2005). 

B. C. Larson, W. Ku, J. Z. Tischler, C.-C. Lee, O. D. Re- 
strepo, A. G. Eguiluz, P. Zschack, and K. D. Finkelstein, 
Phys. Rev. Lett. 99, 026401 (2007). 

M. W. Haverkort, A. Tanaka, L. H. Tjeng, and G. A. 
Sawatzky, Phys. Rev. Lett. 99, 257401 (2007). 
M. van Veenendaal and M. W. Haverkort, Phys. Rev. B 
77, 224107 (2008). 

R. A. Gordon, G. T. Seidler, T. T. Fister, M. W. 
Haverkort, G. A. Sawatzky, A. Tanaka, and T. K. Sham, 
Europhys. Lett. 81, 26004 (2008). 

K. T. Moore and G. van der Laan, Rev. Mod. Phys. 81, 
235 (2009). 

P. Santini, S. Carretta, G. Amoretti, R. Caciuffo, 
N. Magnani, and G. H. Lander, Rev. Mod. Phys. 81, 
807 (2009). 

S. B. Wilkins, R. Caciuffo, C. Detlefs, J. Rebizant, E. Co- 
lineau, F. Wastin, and G. H. Lander, Phys. Rev. B 73, 
060406(R) (2006). 

R. Caciuffo, G. Amoretti, P. Santini, G. H. Lander, 
J. Kulda, and P. d. V. Du Plessis, Phys. Rev. B 59, 
13892 (1999). 

Y. Baer and J. Schones, Solid State Commun. 33, 885 
(1980). 

R. Verbeni, T. Pylkkanen, S. Huotari, L. Simonelli, 
G. Vanko, K. Martel, C. Henriquet, and G. Monaco, J. 
Synchrotron Rad. 16, 469 (2009). 

P. Nozieres and D. Pines, Phys. Rev. 113, 1254 (1959). 
W. Schiilke, H. Nagasawa, S. Mourikis, and P. Lanzki, 
Phys. Rev. B 33, 6744 (1986). 

W. Schiilke, Electron Dynamics by Inelastic X-Ray Scat- 
tering (Oxford University Press, Oxford, UK, 2007). 
G. Kalkowski, G. Kaindl, W. D. Brewer, and W. Krone, 
Phys. Rev. B 35, 2667 (1987). 

M. T. Butterfield, K. T. Moore, G. van der Laan, M. A. 
Wall, and R. G. Haire, Phys. Rev. B 77, 113109 (2008). 
K. T. Moore and G. van der Laan, Ultramicroscopy 107, 
1201 (2007). 

R. D. Cowan, The Theory of Atomic Structure and Spec- 
tra (University of California Press, Berkeley, CA, 1981). 
B. T. Thole and G. van der Laan, Phys. Rev. B 38, 3158 
(1988). 



